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ABSTRACT 


Deep water directional wave spectra, measured by an NDBC 3-meter buoy off 
Montery bay, are transformed to shallow water using a linear refraction model. The 
transformed directional spectra are compared with measured spectra using pressure 
gauge arrays in shallow waters at Marina and Santa Cruz. 

The classical Longuet-Higgins et al. (1963) method of computing directional wave 
spectra and a new exact Fourier coefficients representation method (Grauzinis, 1989) are 
used to compute directional wave spectra. The new method of computing directional 
wave spectra, which represents bi-modal distributions of wave energy exactly matching 
the measured Fourier coefficients to second order, demonstrates improved directional 
resolution over the classical technique. 

This work examines the accuracy and limitations of modeling linear refraction by 
comparing with field observations over complex bathymetry. In general, linear refrac¬ 
tion can give reasonable energy and direction estimates starting with deep water spectra, 
but notable exceptions can occur. The largest prediction error occured at .Marina on 
18 January 1988 for the case of a severe storm. This is presumably due to diffractive and 
non-linear effects of the high waves causing loss of accuracy. The linear refraction 
model is not suitable for handling such problems. 
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I. INTRODUCTION 


A. MOTIVATION OF STUDY 

Wave climatology plays an important role in coastal dynamics. A basic problem is 
how to use deep water wave data to predict shallow water wave characteristics, which 
would allow predictions of waves at any nearshore location using only the measurement 
at one deep ocean location. But there are two associated questions. First . how to char¬ 
acterize ocean wave behavior? Second, how to transform waves from deep water to 
shallow water? Since directional characteristics of ocean waves are important to under¬ 
stand for wave generation and wave refraction behavior, we choose to use the directional 
wave spectral density (DSD) for characterizing the ocean wave behavior and we will use 
a linear refraction model to transform the deep water DSD to shallow water. 

The area of Monterey Bay has been selected for this comparative study of ocean 
waves. This area presents an interesting problem for predicting directional wave spectra 
because of its variable shoreline orientation relative to prevalent eastward waves and the 
presence of the Monterey submarine canyon which is the largest in the western hemi¬ 
sphere and causes substantial perturbations on the nearshore wave field. Figure 1 shows 
the irregrlar bathymetry of the Monterey submarine canyon and three sensors locations 
used in the comparison study. 

A NOAA, 3m buoy is used in deep water to measure the three orthogonal compo¬ 
nents of the wave field heave, pitch, and roll directly. In shallow water, square slope 
arrays, 6m on a side composed of 4 pressure sensors each are used to measure and 
compute the equivalent three orthogonal wave components at Marina and Santa Cruz. 

B. HISTORICAL BACKGROUND 

Historically, Monterey Bay was the first, and presumably one of the most intensively 
studied, locations for wave refraction. Wave refraction diagrams of Monterey Bay were 
first constructed in 1948 by Johnson, O Brien and Isaac (Wiegel, 1964). The first in¬ 
tensive refraction studies were performed for amphibious landing exercises at Fort Ord 
(1954) followed by studies by the Corps of Engineers (1958) for the small craft harbor 
at Santa Cruz. Wave refraction and littoral drift calculations have been performed by 
Dorman (1969) and Amal et al. (1973). All these studies w'ere performed using hand- 
drawm refraction diagrams. More recently, numerical refraction calculations have been 
made for a number of specific marine construction sites (e.g., Thornton (1978, 1980)) 
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and erosion studies (Oradiwe, 1986; McGee, 1986). The present work is an extention 
of the work by Khalid (1989) who used the low’ resolution spectral estimator method of 
Longuet-Higgins, Cartwright and Smith (1963). 

C. OBJECTIVE 

The objective of this study is to investigate W’hether linear refraction can be used to 
transform directional wave spectra from deep to shallow’ water over complex 
bathymetry. Several different ways of computing directional spectral density (DSD) are 
presented. The deep w’ater DSD are transformed to shallow’ water using the refraction 
approach originally proposed by Longuet-Higgins (1957). These transformed spectra are 
compared w’ith DSD measured directly in shallow water for which the data had been 
acquired from pressure gauge arrays located at Marina and Santa Cruz. In Chapter II. 
of this thesis, linear wave theorx’, tw’o methods of computing directional w’ave spectrum, 
and wave spectrum transformation using the theoretical approach will be addressed. 
Chapter III. describes the data acquisition and sensors. Chapter IV. contains the results 
and discussion. Conclusions are presented in Chapter V. 
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II. THEORETICAL BACKGROUND 

A. THE LINEAR WAVE THEORY 
1. Formulation 

Linear surface gravity wave theory’ has been successful in describing many ob¬ 
served wave phenomena. In the formulation of this theory, see for instance Kinsman 
(1965), the basic assumptions can be summarized as 

• wave frequencies, cu, are much greater than the earth's rotation frequency, Q 

• the fluid is homogeneous, or (o^N, the Brunt- Vaisala frequency 

• surface tension effects are negligible, or , where is the frequency of the 

capillary’ waves 

• the fluid motion is irrotational, or V x T = 0, which allows the velocity potential, 
0, to be defined in terms of the water particle velocity vector as K = — 

• the fluid is incompressible, or V • K = 0 

The combination of the assumptions of irrotationality and incompressibility 
leads to the basic equation of linear wave theory, Laplace's equation 

vV=0 (I) 


Analytical solutions to this equation are obtained by applying the following boundary 
conditions: 

• waves are periodic in space and time 

• the bottom boundary is horizontal and impermable, which is prescribed by 

“’U-—^U = 0 (2) 


• the linearized kinematic free surface boundary condition is expressed as 

Cl 


(3) 


and the linearized dynamic boundary condition as 

d<f> , 

where w is the scalar vertical velocity, h is the local water depth , is the instanta¬ 
neous free surface elevation, g is the acceleration of gravity, and t is time. 
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In the linearization process, it is assumed that the wave amplitude a is small when 
compared to the wavelength L and depth, or that and These restrictions 

make the linear wave theory strictly applicable to small amplitude waves of infinitesimal 
steepness. The resulting solution to the wave equation is a single harmonic function. 
Invoking the superposition principle, the instantaneous sea surface elevation can be re¬ 
presented as an infinite sum of sinusoids as expressed by 


oo oe 

Vit) = ^ + £«) ( 5 ) 

/i=l /i=l 

where K„ is the horizontal vector wavenumber, x = is the position vector, tu, is the 
circular frequency and c„ is the initial phase angle. The dispersion relationship is obtained 
from the solution of the wave equation in applying the surface boundary conditions 

(o^ = gk tanh(A:/t) ' (6) 

The spectral component index n is now dropped for convenience. 

Two limiting regions can be considered, deep and shallow water wave approxi¬ 
mations. These simplified solutions correspond to the consideration of asymptotic forms 
for the hyperbolic functions entering the basic solution definitions. These extreme cases 
are generally identified by the ratio ^ , where h is the local water depth and L is the 
wavelengh as computed by the expressions corresponding to each one of the cases. For 
h/L > 1/2 the solution is considered as belonging to the deep water limit and for h/L 
< 1/20 the solution corresponds to the shallow water case. In between the full solution 
must be considered and the solution is known as the intermediate water wave solution. 

The solution of the linear wave equation can be used to describe the wave in¬ 
duced motions. In particular, the relation between w’ave elevation and subsurface 
pressure head p is described by 

open ocean data are invariably collected in the form of time series of processes 
that are inherently random. Spectral analysis techniques are a natural first choice to 
study these data. The previously presented solution to the Laplace's equation gives rise 
to simple spectral relationships between the sea surface elevation spectrum and pressure 
spectra when considering the water column as a constant parameter linear system. Un- 
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der such a system, any two qua. ;ies are connected in the spectral space via a transfer 
function as 


S/cu) = |//(cu)pS,((y) (8) 

where S^(co) and S,(cu) are respectively the output and input energy-density spectra and 
H{(o) is the transfer function. The constant parameter linear system preserves the input 
frequency at the output and modifies the amplitude and phase of each component of the 
input independently. Similar types of relations can be deduced for the cross-spectral 
quantities. For the case of the open ocean wave field and considering the pressure head 
as input and sea surface elevation spectra as output, the applicable transfer function for 
the pressure is 




cosh kz 
cosh k{h + z) 


(9) 


which is the same term as in the brackets of Eq.(7). 

It was noted in the formulation of the linear wave theory that specific phase 
relations exist between the sea surface elevation and pressure. The relative phase be¬ 
tween two quantities in spectral terms is defined as 

Qxy(f^) 

<t>^c) = arc,an-^ (10) 


where 2,,(<u) is the quadrature spectrum and C,(q>) is the co-spectrum of the two 
quantities. 

The spectral transfer function relating pressure head to surface elevation as a 
function of frequency and depth, Eq.(9), can be modified to include differential distance 
and elevation between two sensors of an array to compute energy and cross - spectra 
of surface slope in the x - direction as 



(H) 

(12) 
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where AP, is subsurface pressure difference between two sensors on the x-axis and Ajc is 
differential distance along the same axis. Similarly, pressure power spectra and cross - 
spectra can be computed along y-axis. 

B. THE SPECIFICATION OF A WAVE FIELD 

The complete description of the wave field requires knowledge of frequency and di¬ 
rectional wave spectra. Directional spectra are commonly measured by pitch and roll 
buoy and multi-element arrays, which represent spectra measured at a point. The meth¬ 
ods for computation of the directional wave spectrum are generally classified as 

• .Model fitting methods - based on the parametric representation of the spectrum 
as the classical Longuet-Higgins et al (1963) method. 

• Model independent methods!Davis and Regier,1977), that can be divided into a 
priori and data adaptative methods. 

A priori methods require in advance a definition of criteria to be verified without refer¬ 
ence to the data in use. The data adaptative methods are based in a posteriori assump¬ 
tions about the spectrum , which depart from characteristics of the input data. Any of 
these processes of estimation are conditioned by the traditional constraints inherent in 
time series analysis of statistical reliability and resolution. 

1. Directional Spectrum 

The wavenumber-frequency characteristic of sea waves is contained in the wave 
spectrum as given by the Fourier transform of the correlation function R(r, t) 

S{K, (o) = — i?(r, T)e"'^ • ' dr (13) 

( 2 «) JJ 

This quantity is simply the three dimensional distribution of variance generally known 
as the power spectrum. The knowledge of the wave field at fixed positions (Munk et al., 
1963), permits the computation of the correlation function. The definition of the power 
spectrum as given by Eq.(13) leads by inverse transforming to the correlation function, 
that is 


Ri7, t) = J J5(^, ~^^dKd(o 


(14) 
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which for the assumed stationary wave field is an even function. Consequently, the 
power spectrum is real and symmetric about zero frequency. The correlation-iSanction 
can then be rewritten as a cosine transform of the power spectrum 


R(;, t) = J J SiK, <o) cos{K, (o)dKd(o 


(15) 


The correlation function can now be redefined as 


/?(r,T) = J JsC/l, cu)( cos( cut) cos(Ar«r) + sin(tt)T) sin( AT* (16) 


The power spectrum can then be calculated in terms of cosine and sine trans¬ 
forms of the correlation 


S(k, (o) = jCC(7, w) + iQ(r, co)3e~"^*'dP (17) 

where C(r, <o) and Q(f, oj) are respectively the co- and quadrature spectra of the corre¬ 
lation. Such functions could be calculated if the correlation is known as continuous 
function of the spatial coordinate, although real world observations lead to the col¬ 
lection of data in a small number of positions. Approximate methods of computing di¬ 
rectional spectrum from observations of the <ve field at discrete positions must then 
be considered. 

2. Classical Approach 

Lonquet-Higgins, Cartwright and Smith (1963) following the suggestions of 
Barber (1946), developed a method of computing the directional spectrum from the 
motions of a fioating buoy. The method is based on relationships between the cross 
spectral quantities of the heave, pitch, and roll of the buoy and the Fourier coefficients 
of the directional spectrum. The directional spectrum can be represented as a Fourier 
series 


OO 

S(f,0) * *^ + cos sin nd 

nm\ 


(18) 
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where d is the wave direction and the coefficients a, and are only frequency dependent 
and represent the frequency distribution of energy. Separating the frequency and direc¬ 
tional contributions, as it is usually done, the directional spectrum can be rewritten as 

Sm = Sif)D{f,d) (19) 

where D(f,6) is the spreading function. D{f,6) corresponds to a unity area weighting 
function that redistributes the energy contained in the power spectrum over direction, 
ideally reproducing real world conditions. The three measured quantities obtained with 
a buoy can be directly related with the harmonic representation of the sea surface. The 
heave corresponds to the instantaneous sea surface elevation and the pitch and roll are 
associated with the spatial partial derivatives of the surface elevation. Considering a 
single harmonic of the sea surface representation 

ti{t) = a cos{K • X —(ot) (20) 

the two horizontal cartesian components of the sea slope can be written as 


= 

c>;(r) 

cos 6)fi{t) 

(21) 

= 

^={-iKsxne)r,{t) 

(23) 


Cross spectral quantities can be computed for the instantaneous surface ele¬ 
vation and surface slope, and related with the Fourier coefficients of the directional 
spectrum. These cross spectral quantities can be represented as 


C^{f)^S{f)rD{f,d)de 

•'0 


(23) 


*'0 

Q^JJ)^KS{f)r 

•>0 


D{f,d) cos Odd 


Diffi) sin edB 


(24) 


(25) 
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(26) 


C, ,^(/) = Kh{f) cosHde 

' * **0 

^ (/) = K^SiJ) Sin^ddd (27) 

* * •'0 

(/) = K^S(f) f ^’'z)(/:0) cos 0 sin Bdd (28) 

Considering the directional spectrum in terms of its first five Fourier coefilcients, that 
S(f,d) = + a, cos - 6, sin 6 + ai cos 26 + sin 26 (2‘5) 

, it can be shown that these coefficients are given in terms of the cross spectral quan; 
ties, Eqs.(23-28), by 


— n 

(30) 


(31) 


(32) 


(33) 

bi = Q _ 

(34) 

and equivalently, if a normalization by auto-spectrum is done, the Fourier coefficients 
of the unit area spreading funaion are determined. This approach constitutes the clas¬ 
sical LCS method. The approximate representation of the directional spectrum in terms 

of its Fourier series truncated to its first five coefficients leads to a directional distrib- 


ution of energy as given by a cosine-bell function. Such a representation is of limited 
resolution, and typically has large negative side lobes (Figure 2V 
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To solve the problem of the unrealistic negative side lobes, LCS imposed a 
constraint of positivity on the directional spreading function which resulted in a 
smoothing effect. However this approach has the consequence that the resolution of a 
distribution narrow’er than cos^C —) is not possible. It must be noted that this artifice 
does not solve the real problem that is related to unresolved harmonics of the severely 
truncated Fourier series representation. Also, the LCS approach is not capable of re¬ 
presenting multi-modal distributions of energy that are frequently observed in the real 
ocean, as can be concluded from the truncated Fourier series given by Eq.(29) and 
shown in (Figure 2). 

The method just described can also be applied to measurements using slope ar¬ 
rays, w’hich consist of multi-point arrangements of pressure transducers. Typically the 
sensors in the slope arrays are arranged in a square. The real slope is approximated by 
dividing the differential pressure of a pair of sensors by their separation w’hile assuming 
the sea surface slope as constant between each pair of sensors. The method of com¬ 
puting the directional distribution of energy then proceeds as before for the case of 
measurements with a floating buoy. The sea surface elevation is obtained from the time 
series of pressure by applying the appropriate transfer function. The approximation of 
the surface slope introduces errors in the computation given by (Seymour and Higgins, 
1977) 


dri %in{Kbai) drj 
2 


(35) 


where a, = (cos 6, sin 6) for i= 1,2 and b is the horizontal sensor separation. It can be 
seen that the error of such an approximation is a function of the wavenumber vector and 
the size of the array, decreasing in magnitude as these quantities are reduced. It is then 
natural to employ arrays as small as possible, limited by the resolution and accuracy of 
the pressure sensors. 

3. Exact Fourier Coefllcieiit Representation Method 

The LCS approach provides a way of computing the first five Fourier coeffi¬ 
cients of the series expansion of the directional spreading function. These coefficients are 
obtained via cross- spectral analysis of triorthogonal components of the wave field 
measurements. The number of Fourier coeflicents are limited because of the measure¬ 
ment technique normally used, floating buoy or slope array of four sensors in a square. 
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Figure 2. Typical spreading function of the LCS method: Small directional resol¬ 
ution and unrealistic negative side lobes. 

Such a system can only give information of the wave field through second order terms. 
The coelficicnls determined can be represented by an infinite number of dilferent dis¬ 
tributions whose Fourier expansion is common. Examples of distributions having the 
same 'ouricr coefficients can be seen in Figure 3. It is then natural to choose a phys¬ 
ical!; calisable distribution, which matches the coefficients determined by .observation 
and more exactly simulates the true directional distribution. 

Grauzinis(1989) exploits the matching of a set of Fourier coefficients by dif¬ 
ferent distributions to develop a new method of computing directional spectra. The 
previously defined directional spreading function D{B) can be represented as 


m - [ 1 + cos(*l»,) + bt Sin(*e,)]) (37) 

where a, and fr, are the unitary Ftnirier coefficients. Such coefficients enter the Tinition 
of the polar coefficients c* and as given by 
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(37) 


Q = «* + = <^*^'*^* 

Any symmetric unity area function can be determined solely by its canonical coefficients 
defined as 


rtij^ = I A'(6) cos kddd (38) 

^-ir 

where X{d) is any arbitray unit area function symmetric around 6 = 0. This represen¬ 
tation is simply the cosine transform of an even function. It is seen that a relation can 
be established between the canonical coefficients and the Fourier coefficients. 

Ok +‘h = c*«'*®* = (39) 

where the polar form is now' used for the Fourier coefficients of the spreading function. 
Using linear superposition, any unit area function can be reproduced by a weighted sum 
of symmetric unit area density functions oriented at different directions as given by 

<^*«^*^* ^ (40) 

J 

w’here w, is the weighting factor and Pj is the orientation of each of the components 
considered in the series. Eq.(40) constitutes the generating equation for the polar coef¬ 
ficients where on the right hand side mSj) and Pj are unknowns. The left hand side 
of Eq.(40) is determined from the cross spectral quantities of the wave field typically 
known through second order. A system of equations can be established as 



(41) 


(42) 


which can be seen to be underdetermined. It must be noted that in the system Eq.(4I-42) 
it is assumed that the canonical coefficients are common to both of the distributions of 
energy considered. This was the approach followed by Grauzinis (1989) to obtain 
analytic solutions for the system of generation equations. By the condition of 
orthogonality of the exponentials of different order and the side condition that the sum 
of the weights must be equal to unity, analytic solutions are formed. The analytic 
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DIRECTION 


f jgui c 3. Distribution of energy: The three distribution have a conunon Fourier 
expansion. 

solutions constitute an implicit underdetermined system of four equations. The 
unknowns of this system are the two canonical coefficients, the orientation of the two 
unit area density functions and one of the weighting functions. To solve this system, 
further constraints are needed. Assuming a functional relationship between the 
canonical coefTiceints m, and m, particular solutions are obtained. The original relation¬ 
ship adopted was a power law between the two canonical coefficients 

mj =» m\ (43) 

that for different values of r originated different distributions of energy as analytic sol¬ 
utions. Such distributions are of decreasing sharpness as the magnitude of the exponent 
increases. For the values of the exponent 2, 3 and 4, the resulting distributions in its 
unimodal form are represented in Figure 4 using two different orientations and weights 
also allows this method to represent bimodal distributions of energy. 1 he variable 
sharpness of the model distributions considered will permit the matching of observed 
distributions of considerably more peakedness than the classical LCS method. For a 
canonical coefficient of .9, the haff-power width of the distribution D4 is 12 degrees, that 
is about one-third of the typical resolution of the cosine-bell distribution. 


t4 







Knowledge of Fourier coefficients of order higher than the second introduces 
considerable advantage to this method. Some of the constraints can be dropnfd, such 
as the functional relationship of the canonical coefficients. Moreover, the matching of 
coefficients of higher order increases the resolution of the method. 

The exact matching of the Fourier coefficients through second order of any 
spreading function can be obtained by the method just described. The problem of what 
criteria to use for the selection of the distribution that best matches the real data must 
also be considered. It should be emphasized that there exist an infinity of distributions 
that can be considered to match the same set of Fourier coefficients. The problem of 
selecting the best fitting distribution partly amounts to the matching of the ocean 
peakedness. Cartwright (1963) considered this problem when computing directional 
wave spectra by the LCS method. He concluded that information could be obtained 
from the ratio of the magnitude of the ccefficients C*, Eq.(37). The magnitude of such 
coefficient is given by 

1 Q I = (44) 

Recalling the functional relationship of the canonical coefficients introduced to allow for 
the solution of the generating system, Eq.(43), it is concluded that in the Grauzinis 
method it is more natural to consider the ratio of the natural logarithms of the C* coef¬ 
ficients, that is 


which represents the greatest power-law of the canonical coefficients that can match the 
input Fourier coefficients. Canwright (1963) verified that the ratios obtained from the 
observed data were consistently greater than the values corresponding to the cosine-bell 
function of the LCS method. It was hypothesized that such difference could be due to 
the multimodality of the wind-wave spectrum. Recalling that the ratio of the coefficients 
C» as an indicator was deduced from uni-modal considerations; it does not seem correct 
to make use of it to infer of the adequacy of multi-modal distributions. Grauzinis (1989) 
points out that the natural logarithimic ratio will exceed the power law of canonical co¬ 
efficients increasingly as the beam separation gets larger. Moreover, noise affects the 
coefficients Cj and C, differently. Eqs.(23-28) show that the effects of noise are cancelled 
for the coefficient a,, and consequently the Log Ratio will be decreased. 
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Figure 4. Three model distribution: Uni-modal representation of three din'erent 
model distributions for two values of the canonical coellicient. D2 ,D3 
ind D orrespond respectively to the exponent 2,3 and 4 for the func¬ 
tional i tionship of the canonical coefTicients Eq.(43). I he sharpness 
of the distribution decreases with the magnitude of exponent. For a 
value of the canonical cocfTicient equal to .9 the half power width of D4 
is 12 degrees. 

C. WAVE SPECTRUM TRANSFORMATION 

According to Mehaute and Wang (1982), all the formulations of wave spectrum 
transformation over a shoaling bottom are based on the premise that the wave energy 
associated with a narrow frequency band stays within that band during the transforma¬ 
tion and that it obeys linear superposition. (The shift of energy from frequency to fre¬ 
quency is considered to be the results of nonlinearity.) For each frequency, the energy 
level, which is characterized by the square of a free surface elevation, >f^(0, is treated as 
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an invariant during the transformation. Therefore, the law of tansformation which has 
been established for monochromatic waves is applied as if each value, ?/(f), were identical 
to a periodic wave of the same amplitude and frequency, and that the energy contained 
in each band and direction travels along its corresponding wave ray at group velocity. 

The first rigorous theoretical approach to the evolution of the wave spectrum is due 
to Longuet-Higgins (1957), who demonstrated that the energy spectral density in the 

wave number space remains constant as one follows any wave group along a wave ray 

(Phillips 1966). The demonstration is based on the conservation of energy and the law 
of wave refraction. It is shown that S(/:„/cj) = constant along a wave orthogonal, or 
-^S{ki, ^j) = 0 . Other demonstrations have also been presented by Karlsson (1969), 
Collins (1972), and Krasitskiy (1974) in analogy with geometrical optics and the ray 
particle analogy. 

Following their approach, let us consider a wave spectrum defined in terms of space 

(jc,j'). time(/), and wave number (k„kj) as S{ki, k 2 ,x,y, t). By virtue of the aforemen- 
dS 

tioned principle = 0), then 

cS , cS dx . cS dy ^ 8S dk^ . 6S dk^ ^ 

— —I- - — + -t;- — + -3-h -rr,-7“ = 0 (46) 

0 / ex dt cy dl ok^ dt ok-^ dt 


The first term is generally taken equal to zero under steady- state condition. The second 
and third terms are the convective terms and are due to the variation of energy level as 
a result of the spatial variation. The two last terms are due to shoaling and refraction, 
is the velocity of energy propagation in the x direction and 

^ = (47) 


Similarly: 


dt 


d(o 

dki 


V sin 6 


(48) 


where V: group velocity, k: wave number. Since <u^ = gk tanh kh, oi is a function of k 
and h . The depth, h, depends exclusively on (jf,>'), thus we can write 
(o = coLA,/i,(jc,y)] = (o{k^,y). Then 

du) _ doj dk^ 8co dx , dto dk^ 8o) dy 

dt ck, dt dx dt dk^ dt dy dt 
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in which C «* the phase velocity, C - It is recalled that the ray equations are (Munk 
and Arthur 1951) 


( 56 ) 


ds 


dx a -a 

—T~ = COS 6; “7“ — sin 6 
ds ds 


1 


Y sin64^-cos04^ 

\ cx dy 


in which s = the distance along a ray. Inserting these into Eq.(55) yields simply: 

^^icysm^-o (57) 

i.e.: 

CVS(f,d) = constant (58) 


along a wave ray. Note: 


dS{k„ k,) 
dt 


_V _rf 

Inta ds 


lcvs{r,d)l = o 


(59) 


The result obtained in Eq.(59) is the same as Longuet-Higgins. Eq.(58) can be written 
further as 


— S{f,6rX,y) = constant 


or 


S(co.0) = ^-^So(<u,0o) (61) 

where the subscript refers to deep water quantities. In other words, the ratio between 
shallow w’ater DSD and deep water DSD is a universal function of depth and 
wavenumber where 

(62) 

Figure 5 shows the plots and numerical values of the ratio between deep and shallow 
DSD Eq.(62) at Marina and Santa Cruz versus different frequencies. 
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Figure 5. Truisformation constant along a ray path 
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III. DATA ACQUISITION AND SENSORS 


A. LOCATION 

Monterey Bay was selected as the study area, which poses complex problems due to 
the highly irregular bathymetry of the Monterey submarine canyon. A deepwater wave 
buoy and two shallow water slope arrays were utilized for wave data acquisition. The 
locations of sensors are listed in Table 1, see Figure 1 for bathymetry and location. 


Table 1. LOCATIONS OF SENSORS 


No. 

Type of Sensor 

Location 

Depth (M) 

1. 

NDBC roll; pitch wave 
buoy 

36“48.0N. 

122‘’23.0'59.0"W 

2003 

2. 

Array of pressure gauges 
(4) Marina 

36“42.0'N, 12r48.9'W 

15.0 

3. 

Array of pressure gauges 
(4) Santa Cruz 

36“57.0'N, 122“0.2'W 

13.1 


In comparing coincident data sets, on the average, a time lag of one hour was assumed 
for the waves to travel from the deep water buoy to the shallow water array sites. 

B. DEEP WATER WAVE DATA 

Deep water wave data of the pitch/roll buoy were obtained from the NDBC (Na¬ 
tional Data Buoy Center). The NDBC 3-meter buoy measures heave, pitch, and roll 
movements with respect to a 3-axes magnetometer. The magnetometer axes are aligned 
with fore and after (bow positive), athwartship (starboard positive) and vertical (positive 
upward) axes of a buoy, see Figure 6. 

The analog sensor measures vertical acceleration and displacement, hull pitch and 
roll. The magnetometer bow and starboard output, together with vertical displacement, 
pitch and roll output from analog sensors are sent to the Directional Wave Data Ana¬ 
lyzer (DWDA) after voltage amplification and filtering. 

The hull azimuth angle is calculated from bow and starboard magnetic components 
with necessary corrections made for pitch, roll, and hull magnetic fields. The magnetic 
azimuth is converted into true azimuth using the magnetic variation for the site. The 
true azimuth is then used to calculate East-West and North-South components of hull 
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Figure 6. Floating Buoy 


slope. One hertz time record of hull displacement and slope are stored in DWDA for 
spectral analysis and transmission to shore station through GOES (Steele et al, 1985). 

Data provided on magnetic tape included environmental parameters, wave spectra, 
directional wave parameters, and co and quad-spectra. Directional wave parameter data 
were utilized to determine directional wave spectrum. 

C. SHALi-OW WATER WAVE DATA 

The calibrated data for the presst s gauge arrays were obtained from the Scripps 
Institute of Oceanography. Array consists of four pressure gauges located in a six-meter 
square area which measures pressure head in centimeters of water, See Figure 7. Each 
sensor is located about 50 to 100 cm above the bottom. Pressure power spectrum is 
calculated which in turn is transformed into surface elevation fi and surface slopes tj, and 
spectra by applying linear wave theory transfer function as described in Chapter II. 
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Figure 7. Slope Pressure Array 


For the computation of the cross-spectra, the following values were utilized: 

• sampling interval = At “ 2 sec 

• total data points ■» 2048 

• data points for one sample » 64 

• nyquist frequency * 0.25 Hz 

• frequency bandwidth - 4/^ - 0.0078125 Hz 

The spectra were cut at a frequency of 0.2 Hz in the higher frequency range due to 
overestimation by the transfer function. The wave‘period for high quality estimate is 
approximately 5-20 seconds. 




D WAVE DATA SELECTION 

Selected wave cases from January 1988 were used in the analysis. During the winter 
season, weather conditions were generally favorable for creating higher energy waves. 
The selected days are 3 January 0800 PST, 18 January at 0200 PST and 29 January^ at 
2100 PST. These periods provide two cases of narrow band energy spectra at different 
frequencies and directions and a third case of a major storm event. The California coast 
w'as hit by an unusually severe storm during 16-18 January 1988 which resulted in high 
waves, strong winds, and sea level well in excess of predicted values (Cayan, et al., 1988). 
Widespread damage occured along the southern California coast. This storm developed 
about 300 XM west of San Francisco on 16 January 1988 traveled south-east at about 
33 knots. The storm center passed off Monterey Bay on the morning of 17 January 1988 
continuing on its SE track and making landfall at Avila Beach in Central California. 
The energy spectra for these cases are depicted in Figure 8, 9, and 10 comparing the deep 
water buoy and shallow water pressure gauge measurements at Marina and Santa Cruz. 
Deep water wave characterics of selected cases are given in Table 2. 


Table 2. DEEP WATER WAVE CHARACTERISTICS OF SELECTED CASES 


Characteristics 

3 Jan 88 

18 Jan 88 

29 Jan 88 


1.87 m 

5.12 m 

2.65 m 

Freq,„, 

0.13 Hz 

0.07 Hz 

0.09 Hz 


0.18 Hz 

0.10 Hz 

0.12 Hz 

Direction„,„ 

o 

O 

oo 

300° 

275° 


Reviewing Table 2 and other published data, deep water waves generally travel 
eastward during the month of January with slight variability of about 20° on eigher side 
of true East. This direction is dominant in the energetic middle frequency range of 0.08 
to 0.13 Hz. At other frequencies, the direction tends to be variable. 

The significant wave heights, calculated as H, = AyJa ^, where is variance calcu¬ 
lated as the area under the energy density spectra, are given in Table 3. 

E. TRANSFORMATION OF DEEP WAVE SPECTRUM 

The directional spectrum density at wavenumber space is constant along ray path. 
Transforming from wavenumber space to frequency and direction space, the shallow 
water DSD is related to the deep water DSD by a universal function of depth, h, and 
wavenumber, k, or frequency,/, Eq.(62), however, the relationships betw’een the shallow 
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Figure 8. Measured spectral densities on 3 January 1988, 0800 PST 
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Figure 10. Measured spectral density on 29 January 1988, 2100 PST 
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Table 3. SIGNIFICANT WAVE HEIGHT (UNIT: METER) 



Buoy 

Marina 

Santa Cruz 

3 Jan 1988 

1.87 .M 

0.8 .M 

0.9 M 

18 Jan 1988 

5.12 .M 

5.15 M 

2.4 M 

29 Jan 1988 

2.65 .M 

1.77 .M 

0.94 .M 


water angle and deep water angle are still to be found. The Dobson linear refractir 
model was used for this purpose. The reverse projection method was employed. That is, 
the linear refraction model program was run backward from shallow water location to 
deep water, beginning at the shallow water locations at Marina and Santa Cruz. Figure 
12 - 15. show rays at the particular frequencies of 0.03, 0.06, 0 ). 0.13 Hz from .arina 
and Santa Cruz propagated offshore in increments of 0.5 degrees over the range of all 
possible incoming wave angles. Once the de. o water is reached, the rays are stopped 
and deep water wave directions are recorded. As can be seen the low frequency waves 
are more refracted by the bathymetry than high frequency waves. 

An example of the shallow w’ater angle versus deep water angle at Marina for 0.06 
Hz is show'n in Figure 15. The doited line is the least square first order fit, where a is 
the slope and b is the intersection at shallow water angle axis when deep water angle 
becomes zero. 

The slope and intercept of least squa. inear fits at various frequencies and lo¬ 
cations are given in Table 4 and plots given a .Appendix A. As frequency increases, the 
slope increases to one and the offset decreases toward zero. Since slope less than one 
amounts to shrinking the deep water angular range, it indicates that smaller slope will 
concentrate wave energy into a smaller direction range and decrease the spectral density. 
The offset is equivalent to rotating the deep water DSD to a different direction. 
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Figure 12. Ray traces from Marina to deep water 
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Figure 13. Ray traces from Santa Cruz 
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Shallow Water Angle vs Deep Water Angle at Marina, 0.06 Hz 



Figure 15. Example of the shallow water angle versus deep water angle 






Table 4. EFFECTS OF LINEAR REFRACTION _ 

Marina Santa Cruz 


Frequency 

slope 

olTset 

slope 

offset 

0.03 Hz 

0.1164 

242.0 

0.3290 

126.9 

0.06 Hz 

0.3020 

189.0 

0.3974 

116.6 

0.09 Hz 

0.4165 

157.9 

0.4998 

98.0 

0.13 Hz 

0.6938 

82.8 

0.6824 

62.1 
























IV. RESULTS AND DISCUSSION 


A. COMPARISON OF LCS AND EFC MEASURED DIRECTIONAL SPECTRAL 
DENSITIES 

A flow chart of the calculation of measurements is shown in Figure 16. The 
TRIDEN program calculates the EFC spectrum and compares it with the spectrum 
measured in shallow water. The directional distribution of measured waves in polar plot 
form are shown in Figures 17,18, and 19 at the three locations. Buoy, Marina, and Santa 
Cruz on 3 January 1988. These axis are aligned with the earth directions. The deep 
water wave approach from northwest or west. Four frequencies, 0.03,0.06, 0.09 and 0.13 
Hz are examined. At Marina, waves come mostly from west, and at Santa Cruz, from 
southwest. The plots show the energy is more widely distributed in deep w'ater and that 
the spectrum narrow's in shallow water, except in the low'est frequency. 

The energy density distribution as a function of frequency and direction are shown 
in Figures 20 - 31, for 3 January 1988 at deep water Buoy, Marina, and Santa Cruz. 
Both of the LCS method and EFC method indicate the same directionality of measured 
data, but EFC method show’s a much sharper peak than LCS method. Also LCS results 
sometimes show negative side lobes, which do not exist in the EFC results. 

The wave energy spectra at frequency 0.13 Hz, are shown in Figures 32, 33, and 34. 
Again , the EFC method results in much higher resolution. For this reason, only the 
EFC are discussed for other results. The other two days, 18 and 29 January 1988, 3-D 
and 2-D also show similar results. All graphs are given in Appendix B. The deep water 
buoy and shallow water array (Marina and Santa Cruz) energies are compared in Table 
5, 6 for the different methods. As can be seen the LCS and EFC methods give the same 
total energies under the spectra. 


Table 5. DEEP WATER BUOY ENERGY (AP) 


DATE 

EFC 

LCS 

D2 

COMPOSITE 

3 Jan 88, 0800 

0.22 

0.22 

0.22 

18 Jan 88, 0200 

1.642 

1.641 

1.642 

21 Jan 88, 2100 

0.443 

0.442 

0.443 
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Figure 16. Tlte flotr chart of measurement calculation 


Table 6. SHALLOW WATER ARRAY ENERGY (AP) 


DATA 

EFC 

LCS 

D2 

COMPOSITE 

Marina 

Santa 

Cruz 

Marina 

mmrm 

Marina 

Santa 

Cruz 

3 Jan 88.080U 

0.042 

0.054 

0.042 

0.054 

0.042 

0.054 

18 Jan 88, 
0200 

1.659 

0.363 

1.656 

0.362 

1.656 

0.362 

21 Jan 88. 
2100 

0.196 

0.057 

0.196 

0.057 

0.196 

0.057 


B. REFRACTION MODEL PREDICTIONS 

The deep water spectrum is transformed to shallow water using Eq.(62). 1 he shal- 
low water spectrum is simply equal to deep water spectrum multiplied ’^y shoaling. 
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BUOY 


Figure 17. The energy direction of deep water on 3 Jan 1988. 
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SANTA CRUZ 


Figure 19. The energy direction of shallow water at Santa Cruz on 3 Jan 1988. 
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Figure 20. The LCS directipnel spectral density in 3-D plot: On 3 Jan 1988 at 
deep water Buoy. 
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OlRECnONAL SPECTRA 
Deep Water Buoy 
Jan 03, 1960 0800 POT 



Figure 21. The EFC directional spectral density in 3-D plot: On 3 Jan 1988 at 
deep water Buoy. 
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Figure 22. The LCS directional spectral density' in contour plot: On 3 Jan 1988 

at deep water Buoy. 
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Figure 23. The EFC directional spectral density in contour plot: On 3 Jan 1988 
at deep water Buoy. 
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Figure 24 . The LCS directional spectral density in 3-D plot: On 3 Jan 198S at 
Marina. 
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Figure 25. 


The EFC directional spectral density in 3-D plot; On 3 Jan 1988 at 
Marina. 
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Figure 27. The EFC directional spectral density in contour plot: On 3 Jan 1988 
at Marina. 
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Figure 28. The LCS directional spectral density in 3>D plot: On 3 Jan 1988 at 

% 

Santa Cruz. 








Figure 29. 


The EFC directio|ial spectral density in 3-D plot: On 3 Jan 1988 at 
Santa Cruz. 
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Figure 31. The EFC directional spectral density in contour plot: On 3 Jan 1988 
at Santa Cruz. 
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Figure 32. The LCS & EFC directional spectral densities in 2-D plot: On 3 Jan 
1988 at deep water Buoy. 
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Figure 33. The LCS &. EFC directional spectral densities in 2-D plot: On 3 Jan 
1988 at Marina. 
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Figure 34. The LCS St EFC (directionai spectral densities in 2-D plot; On 3 Jan 
1988 at Santa Cruz. 
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Table 7. COMPARISON OF TRANSFORMED AND MEASURED DATA DIF¬ 


FERENCE IN DEGREE 


Location 

Marina 

Santa Cruz 

Frequency 

0.06 

Hz 

0.09 

Hz 

0.13 

Hz 

0.06 

Hz 

0.09 

Hz 

0.13 

Hz 

03 Jan 19S8 

25“ 

0“ 

0“ 

15“ 

0“ 

30“ 

18 Jan 1988 

30“ 

95“ 

50“ 

5“ 

15“ 

25“ 

29 Jan 1988 

5“ 

0“ 

5“ 

40“ 

5“ 

5“ 


Table 8. MEASURED AND TRANFORMED ENERGY SPECTRUM DATA AT 


MARINA (APIHZ) 


FREQ. 

0.06 Hz 

0.09 Hz 

0.13 Hz 1 

.MEAS¬ 

URED 

TRANS¬ 

FORMED 

MEAS¬ 

URED 

TRANS¬ 

FORMED 

MEAS¬ 

URED 

TRANS¬ 

FORMED 

3 Jan 88. 0800 

0.026 

0.123 

0.062 

0.078 

0.401 

0.449 

18 Jan 88. 0200 

12.8 

3.85 

1.54 

2.46 

4.11 

2.90 

29 Jan 88. 2100 

0.103 

0.103 

1.28 

1.28 

0.180 

0.205 


Table 9. MEASURED AND TRANFORMED ENERGY SPECTRUM DATA AT 
_SANTA CRUZ {APIHZ) 


FREQ. 

0.06 Hz 

0.09 Hz 

0.13 Hz 1 

MEAS¬ 

URED 

TRANS- 

FOR.MED 

MEAS¬ 

URED 

TRANS¬ 

FORMED 

MEAS¬ 

URED 

TRANS¬ 

FORMED 

3 Jan 88. 0800 

0.072 

0.187 

0.045 

0.II8 

0.205 

0.277 

18 Jan 88. 0200 

2.87 

3.25 

1.03 

2.00 

0.411 

1.54 

29 Jan 88. 2100 

0.062 

0.16 

0.616 

1.03 

0.168 

0.205 


{-p'), and dispersion, {-j^) and the energy density is replotted based on the angle con¬ 
version between deep and shallow water, determined by the linear wave refraction model. 

The Marina location is shadowed by Point Finds in the southwest and Point Santa 
Cruz in the northwest. High frequency waves, O.ll Hz and above, are limited to di¬ 
rections from 235®-310'’ at the array site, whereas low frequency waves may arrive from 
even southerly directions after considerable refraction (refer to Figures 11, 12). 
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The Santa Cruz site poses more compif^x problems than Marina because of its lo¬ 
cation. The waves approaching from the westerly direction ■ ke almost a 90° turn to 
arrive at the shoreline. The low frequency waves travelling ; rth and north-eastward 
may arrive at this location. The ray patterns show (Figures 13, 14) scattered and cross¬ 
ing rays (caustics) from the west or higher approach angle at low frequency. The ray 
patterns become more regular with increasing frequency. 

C. COMPARISON OF DEEP, SHALLOW AND TRANSFORMED 
DIRECTIONAL SPECTRAL DENSITIES 

Transformed and measured shallow wave directional spectra for Marina are com¬ 
pared in Figure 35 - 37. On 3 January 1988 for 0.06 Hz, the result shows abou» a 25 
degree angle difference between transformed and measured DSD, and the transformed 
spectral energy is overestimated. For 0.09 Hz, the transformed and measured shallow 
waves have the same directional spectral distribution, but the measured data has greater 
energy density than the transformed. On the same day at 0.13 Hz, the result shows good 
prediction in directionality. Both frequencies indicate the transformed DSD have nar¬ 
rower distributions compared with deep water DSD, because the shinking of the angle 
distribution concentrates the energy in shallow water during transformation. The re¬ 
fraction effects are greatest at lower frequency. On 18 January 1988 for 0.09 Hz at 
.Marina, the results (Figure 38) show a poor prediction, the directionality having a 95 
degree difference, which is at low relative energy within the spectrum and appears to be 
a spurious result. On 29 January 1988, the results show a 5 degree difference or less at 
all frequencies. 

Results of the Santa Cruz location are dip’ 'ed in Figure 40 - - The results ar: 
different from Marina. On 18 January 1988 at i*9 Hz, the predictic as a 15 degree 
difference, which is better than that at Marina, m 3 January at 0.09 Hz and 29 January' 
at 0.13 Hz both provide satisfactory results. Table 7 summerize the three days 3, 18, and 
29 January 1988 at different frequencies and locations comparing transformed and 
measured directionality angle differences. 

The variance within half power band width of the directional spectra at Marina and 
Santa Cruz on all these days are given in Tables 8 and 9. By examing Tables 7 - 9 and 
focusing on the energy peak frequencies for the three days, 0.13 Hz on 3 January 1988, 
0.06 Hz on 18 January 1988 and 0.09 Hz on 29 January 1988, we can compare the en¬ 
ergy and directional difference between measured and transformed data at .Marina a: J 
Santa Cruz. At Marina, the directions are the same on 3 and 29 January 1988, but has 
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S(f.tltcia). M**2/ll7yRad 


Deep, Shallow, and Transfonned DSD at Marina on 1/3,0.09 Hz 



Figure 35. The deep, shallow, and transformed DSD; On 3 Jan 1988 at Marina, 
0.09 IIz. 
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S(f,ilieia), 


Deep, Shallow, and Transfomed DSD at Marins on 1/3,0.13 Hz 



Direction, Deg 


Figure 36. The deep, shallow, and (ransfomied DSD: On 3 Jan 1988 at Marina, 
0.13 Hz. 
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I//Kn(l 


Deep, Shallow, and Transformed DSD at Marina on 1/3,0.06 Hz 






Deep, Shallow, and Transfonned DSD at Marina on 1/18,0.09 Hz 



Figure 38. The deep, shallow, end transformed DSD; On 18 Jan 1988 at Marina, 
0.09 Hz. ' . 
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S(f.llicta), M**2/H7yRad 


Deep, Shallow, and Transformed DSD at Marina on 1/29,0.13 Hz 



Figure 39. The deep, shallow,.and transformed DSD: On 29 Jan 1988 at Marina, 
0.13 Hz. 


61 











a 30 angle degree difference on 18 January 1988, on the same day, th ' easur crgy 
is four times larger than predicted, but on 3 and 29 Januar>’ 1988 e. s arc . same 
for between measured and predicted. At Santa Cruz, the directional crence is almost 
the same on 18 and 29 January 1988, but has a 30 angle degree difference on 3 January 
1988. From an energy point of view, 3 and 18 January 1988 show the same energy, but 
on 29 Januarv’ 1988, the predicted energy is almost two times the measured energy. 

Summerizing, the energies at the peak frequency were similar except at Marina on 
the 18th (x4) and Santa Cruz on the 29th (x2). The angles at the peak frequencies were 
within 5 degrees except at Santa Cruz on the 3rd (30°) and Marina on the 18th (30°). 
There is no obvious correlation with error on deep water wave directions (~275° on the 
3rd, 300° on the 18th and 245° or. :lie 29th) 
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Deep, Shallow, and Transformed DSD at Sanu Ciuz on 1/18,0.09 Hz 



Figure 40. The deep, shallow, and transformed DSD: On 18 Jan 1988 at Santa 
Cruz, 0.09 Hz. 
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S(f,tlieia), M»*2/ll7yRad 


Deep. Shallow, and Transfonned DSD at Santa Ciuz on 1/3,0.09 Hz 



Figure 41. The deep, shallow^ and transformed DSD: On 3 Jan 1988 at Santa 
Cruz, 0.09 Hz. 
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V. CONCLUSION 


Deep water directional wave spectra were calculated from data acquired from the 
NDBC buoy and transformed over the highly complex bathymetry of Monterey Bay to 
shallow water using a linear wave refraction model. The transformed spectra were com¬ 
pared with measured shallow wave spectra at .Marina and Santa Cruz. 

Two methods of calculating the directional spectrum were employed, the LCS and 
EFC techniques. The results show that the two techniques give the same peak densities 
and energies, but that the EFC method has much superior resolution. 

Tl' - comparison of measured and transformed spectra show that higher frequency 
wave c, .nditions exhibit more accurate results, which are less refracted, and with rela¬ 
tively larger errors occurring at lower frequencies. In general, linear refraction can give 
reasonable energy and direction estimates starting with deep water spectra, but notable 
exceptions can occur. The largest prediction errors occured at .Marina on 18 January 
1988 for the case of a severe storm. This is presumably due to diffractive and non-linear 
effects of the high waves causing loss of accuracy. The linear refraction model is not 
suitable for handling such problems. Further modeling efforts for Monterey Bay includ¬ 
ing diffractive and non-linear effects are needed. 
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Shallow Water Angle, 


APPENDIX A. THE SHALLOW WATER ANGLE VERSUS DEEP 

WATER ANGLE 


Shallow Water Angle vs Deep Water Angle at Marina, 0.09 Hz 



Deep Water Angle,' Deg 


Figure 43. The shallow water angle versus deep water angle at Marina 0.09 llz. 
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Slialiow Water Angie, 



Figure 44 . The shallow water angle versus deep water angle at Marina 0.13 Ilz. 
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Shallow Water Angle. 


Shallow Water Angle vs Deep Water Angle at Santa Cruz, 0.06 Hz 



Figure 45. The shallow water angle versus deep water angle at Santa Cruz 0.06 I Iz. 
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Shallow Water Angle, 


Shallow Water Angle vs Deep Water Angle at Santa Cruz, 0.0° Hz 
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Figure 46 . The shallow water angle versus deep water angle at Santa Cruz 0,09 Hz. 









Shallow Water Angle, 



Figure 47. The shallow water angle versus deep water angle at Santa Cruz 0.13 Hz. 







APPENDIX B. THE DEEP, Si ALLOW, AND TRANSFORMED DSD 


Deep, Shallow, and Transfonned DSD at Marina on 1/18,0.06 Hz 



r'':;ure 48. The deep, shallow, and transformed DSD: On 18 Jan 1988 at Marina, 
0.06 Hz. 


72 








Deep, Shallow, and Transformed DSD at Marina on 1/18,0.13 Hz 



Figure 49. The deep, shallow, and transformed DSD: On 18 Jan 1988 at Marina, 
0.13 Hz. 
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Figure 51. The deep, shallow, and transformed DSD: On 29 Jan 1988 at Marina, 
0.09 Hz. 
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S(f.thcta), M**2/Hz/Rad 
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Figure 53. The deep, shallow, and transformed DSD: On 3 Jan 1988 at Santa 
Cruz, 0.13 Hz. 
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Deep, Shallow, and Transformed DSD at Santa Cruz on 1/18,0.06 Hz 



Figure 54. The deep, shallow, and transformed DSD: On 18 Jan 1988 at Santa 
Cruz, 0.06 Hz. 
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Deep, Shallow, and Transformed DSD at Santa Cruz on 1/18,0.13 Hz 



Direction, Deg 


Figure 55. The deep, shallow, and transformed DSD: On 18 Jan 1988 at Santa 
Cruz, 0.13 Hz. 
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S(f,tljeta), M**2/Hz/Rad 


Deep, Shallow, and Transfoimed DSD at Santa Cruz on 1/29,0.06 Hz 



Figure 56. The deep, shallow, and transformed DSD: On 29 Jan 1988 at Santa 
CTuz, 0.06 Hz. 
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Direction, Deg 

Figure 57. Tlie deep, shallow, and transformed DSD: On 29 Jan 1988 at Santa 
Cruz, 0.09 I Iz. 
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